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We obtain time dependent q-Gaussian wave-packet solutions to a non linear Schrodinger equation 
recently advanced by Nobre, Rego-Montero and Tsallis (NRT) [Phys. Rev. Lett. 106 (2011) 10601]. 
The NRT non-linear equation admits plane wave-like solutions (g-plane waves) compatible with the 
celebrated de Broglie relations connecting wave number and frequency, respectively, with energy 
and momentum. The NRT equation, inspired in the g-generalized thermostatistical formalism, is 
characterized by a parameter g, and in the limit q — > 1 reduces to the standard, linear Schrodinger 
equation. The g-Gaussian solutions to the NRT equation investigated here admit as a particular 
instance the previously known g-plane wave solutions. The present work thus extends the range of 
possible processes yielded by the NRT dynamics that admit an analytical, exact treatment. In the 
q — > 1 limit the g-Gaussian solutions correspond to the Gaussian wave packet solutions to the free 
particle linear Schrodinger equation. In the present work we also show that there are other families 
of nonlinear Schrodinger-like equations, besides the NRT one, exhibiting a dynamics compatible 
with the de Broglie relations. Remarkably, however, the existence of time dependent Gaussian-like 
wave packet solutions is a unique feature of the NRT equation not shared by the aforementioned, 
more general, families of nonlinear evolution equations. 



I. INTRODUCTION 



A nonlinear Schrodinger equation has been recently advanced by Nobre, Rego-Monteiro and Tsallis jUE]. The NRT 
proposal constitutes an intriguing contribution to a line of enquiry that has been the focus of continuous research 
activity for several years: the exploration of non-linear versions of some of the fundamental equations of physics [3111]. 
The NRT equation is inspired in the thermostatistical formalism based upon the Tsallis S q non-additive, power- 
law entropic functional, whose applications to the study of diverse physical system and processes have attracted 
considerable attention in recent years (see, for instance, [SH5] and references therein). In particular, the S q entropy 
constitutes a useful tool for the analysis of diverse problems in quantum physics [TfJlfTT] . 

The NRT nonlinear Schrodinger equation governing the field &(x,t) ("wave function") corresponding to a particle 
of mass m reads [U [2] , 
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where the scaling constant $q guarantees an appropriate physical normalization for the different terms appearing in 
the equation, i is the imaginary constant, ft is Planck's constant, and q is a real parameter formally associated with 
Tsallis' entropic index in non-extensive thermostatistics [6]. It was shown in [T] that the wave equation ([I]) admits 
time dependent solutions having the "g-plane wave" form, 



$(x,t) = $ [1 + (1 - q)i(kx - wt)] 



(2) 



with k and w real parameters having, respectively, dimensions of inverse length and inverse time (that is, k can be 
regarded as a wave number and w as a frequency). In the limit q — > 1, the g-plane waves ^ reduce to the plane wave 
solutions $0 ex P( — i{kx — wt)) of the standard, linear Schrodinger equation describing a free particle of mass m. 

The g-plane wave solutions ([2| propagate at a constant velocity c — w/k without changing shape, thus exhibiting 
a soliton-like behaviour. Moreover, and in contrast to the q — > 1 case yielding standard plane waves, the solutions 
corresponding to q ^ 1 don't have a spatially constant modulus. In fact (defining ip = Q/$>o) we have, 



\^(x,t)\ 2 = [l + il-qfikx-wt) 2 )] 1 - 



(3) 



which corresponds, for 1 < q < 3 to a normalizable g-Gaussian centered at x — wt/k. Therefore, in this case the 
g-plane wave solution describes a phenomenon characterized by a certain degree of spatial localization. A field- 
theoretical approach to the NRT equation was developed in [5] , where it was shown that this equation can be derived 
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from a variational principle. The nonlinear NRT equation is formally related to the nonlinear Fokker-Planck equation 
(NLFP) with a diffusion term depending on a power of the density. These kind of evolution equations, and their 
relations with the nonextensive thermostatitical formalism, have been the focus of an intensive research recently |18l - 
|2"T] . In spite of the formal resemblance between the NRT Schrddinger equation and the nonlinear Fokker-Planck, there 
are profound differences between these two types of equations. For instance, the nonlinear Fokker-Planck equation 
does not admit g-plane wave solutions of the form ([2]), that propagate without changing their shape. 

A property of the solutions ^ that was highlighted by NRT [T] is that they are consistent with the celebrated de 
Broglie relations [2"5] . 

E = hw, 

P = hk, (4) 

connecting, respectively, energy with frequency and momentum with wave number. Indeed, it can be verified that 
the g-plane wave |2| satisfies the equation ([!]) if and only if the parameters w and k comply with the relation, 

hk 2 

w = i (5) 
2m' w 

which, combined with Q, lead to the standard relation between linear momentum and kinetic energy, 

E = h < 6 > 

This suggests that it is conceivable that the g-plane wave ^ represents a particle of mass m with kinetic energy hw 
and momentum hk [1]. 

Wave packets (in particular Gaussian wave packets) are of paramount importance in quantum mechanics, both from 
the conceptual and practical points of view [28, 29 . Wave packets also played a distinguished role in the historical 
development of quantum physics 30J . It is interesting to explore the existence of time dependent wave packet solutions 
of the NRT nonlinear Schrddinger equation. The original presentation of the NRT equation [1] was strongly focused 
upon the g-plane wave solutions. However, as a first step towards elucidating the meaning of the NRT equation, it is 
imperative to explore more general solutions. The aim of the present effort is to investigate a family of exact analytical 
time dependent solutions to the NRT equation exhibiting the form of g-Gaussian wave packets and corresponding, in 
the limit q — > 1, to the celebrated Gaussian wave packets solutions of the linear Schrddinger equation. 



II. TIME DEPENDENT WAVE PACKET SOLUTIONS 



In this work we are going to investigate solutions to the NRT equation based upon the g-Gaussian wave packet 
ansatz, 



ip(x,t) = 



$ 



= [l-(l-q)(a(t)x 2 + b(t)x + c(t))]^ , 



(7) 



where a, b, and c are appropriate (complex) time dependent coefficients. Notice that ip depends on time only through 
these three parameters. Inserting the ansatz Q into the left and the right hand sides of the NRT equation ([T]) one 
obtains, 



ih^- = ~ih(a(t)x 2 + b{t)x + c{t))ip q , 



and, 



1 d 2 ip 2 -" 
2 — g dx 2 



2(3 - g)a(t) V - 2(3 - q)a{i)b{t)x + 2a(t) - 2(1 - q)a{t)c(t) - b(t)' 



(8) 



(9) 



Combining now equations ([lj, ([8]), and ^ one sees that the ansatz ([7 1 constitutes a solution of the NRT equation 
provided that the coefficients a, b, and c comply with the set of coupled ordinary differential equations, 



ia(t) = —(3-q)a(t) 2 
m 

ib{t) = ~(3-q)a(t)b(t) 
m 
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(12) 



3 



The above set of differential equations admits the general solution, 
1 

a(t) 
b(t) 
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m 

(3 



a 



c(t) 



(3 — q)iht 



(3 — q)ih 



t + a 



1-3 
"3-, 



3-9 



1-9 



/3 2 /(3-«)*n \ 3 -« i 

— — i + a +7 

4 V y ' 1-g 



where a, /3, and 7 are integrations constants determined by the initial conditions a(0), 6(0), and c(0), 
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A. Limit Case q — > 1 

Let us now briefly consider the limit q — > 1 of the evolving g-Gaussian wave packet. In this case the time dependent 
parameters a(t), b(t), and c(t) are, 
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with a, /3 and 7 integration constants. The general solution can be written as 



ip = lim \l- (l-q)(a(t)x 2 + b(t)x + c(t))] 1 - q 

q-¥l 

( 1, (2iU \ {/3/2 + x) 2 N 



(22) 



Taking now exp(47) = ^exp(— a) and defining fc = i/3/a, and tan2(9 = (22 1 can be cast under the guise, 



2a/ir 



1/4 



exp 



^— i (6 + hk 2 ) t/2m))^J exp(ikox) 
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(23) 



recovering the well known Gaussian wave packet solution of the standard linear Schrodinger equation. 
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B. Illustrative Example of Wave Packet Evolution for q — 2. 



The evolution of the time dependent g-Gaussian solution is illustrated in Figure 1, where the square modulus \ip\ 2 of 
an initially localized solution is depicted against the (nondimensional) time variable i — f^f) t and spatial coordinate 
x = ^t=£, for q — 2 (here ao = |a(0)| stands for the modulus of the intial value of the parameter a). It is interesting 

that in the nonlinear (q = 2) case, as time progresses, t) | 2 develops two peaks that depart from each other. This 
behaviour exhibits some qualitative similarity with the evolution of an initially localized particle in the tight binding 
model [5T] . 




FIG. 1: Behaviour of |^(a;,£)| 2 determined by the NRT nonlinear Schrodinger equation with q = 2 and initial conditions given 
by a = 1, j3 = 1, and 7 = 1. All depiced quantities are dimensionless. 



C. The Case q = 3 and "Frozen" Solutions 



When q = 3 the set of coupled differential goberning the evolution of the parameters a, b, and c, admit the solution, 



a = a c 
b = b r 



bl - 2a c 



+ c\ exp [2i— a c t , (24) 



where a c , 6 C , and c±, are time independent constants. The concomitant solution to the NRT equation reads, 



ip(x,t) 



1 

7^ 



a c x + b c x + 



4dr. 



ci exp I 2i — a c t 
m 



(25) 



It is interesting that the solution (25) is periodic (with period, in the dimensionless time variable t — —a c t, equal 
to 7r) even though the particle is not subjected to an external confining potential. Indeed, this solution describes a 
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quasi-stationary scenario where the shape of \ip(x,t)\ 2 "pulsates" with the abovementioned period. In the extreme 
case given by c\ = 0, the amplitude of the "pulsations" vanishes and we obtain the stationary solution, 



1 

71 



a c x + b c x + 



4a r „ 



(26) 



The norm N = J \ip\ 2 dx of this last solution is finite provided that a c ^ and b c /a c is not a real number. 

It is interesting that there are "frozen" solutions like ( [26] ) also for instances of the NRT equation characterized by 
other values of q. Indeed, it can be verified that the NRT equation admits stationary solutions of the form 



V>(x) = [bx + ic] 2 -i , (27) 
with b 7^ and c ^ real constants. The square modulus |?/;| 2 of the "frozen" solutions (27 1 has a (/-Gaussian profile, 



\iP{x)\ 2 = [b 2 x 2 + c 2 ] 2 ~* (28) 

and a finite norm for 2 < q < 4. Notice that the solutions ([27]) can be cast under the guise of a ^-exponential, but 
with a q- value given by q = q — 1, which is different from the q- value associated with the concomitant nonlinear NRT 
evolution equation. Therefore, strictly speaking, the solutions (27) are not comprised within those corresponding to 
the ansatz ([7]). 



D. g-Plane Waves and Related Particular Solutions 



The case a = leads to a family of particular solution of the set of equations (10), given by, 
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where b c and cq are complex constants. A simple but important case is obtained when the constant b c is purely 
imaginary and cq = 0, 



= 0, 



—ik, k € R, 

hk 2 



(30) 



which, after setting w = hk 2 /2m, is clearly seen to correspond to the q-plane wave solutions considered in pQ. 
Therefore, the exact solutions to the NRT equation studied in [1] constitute a particular case of the time dependent 
g-Gaussian wave packet investigated here. 



Another interesting particular instance of the solutions associated with ( 29 1 corresponds to the case where the 



constant b c is a real number and cq — —i^^to, with to > a real constant with dimensions of time, leading to, 



ih 



l + (l-g)( -b c x+^-b 2 c (t + t ) 



ip(x,t) 

The squared modulus of the above soution has a g-Gaussian form, 



(31) 



\TP(X, t)\ 2 = 







1*1 





2m(l - (1 - q)b c x) 
(l-q)hb 2 c (t + t ) 
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(32) 



G 



leading to a finite norm, N = J \4>\ 2 dx, for t > —to and 1 < q < 3 (see next subsection). The solution (31) exhibits 
a finite-time singularity in the past: indeed, its norm diverges at the finite time t = —to- In the limit q — > 1 the 
solutions pll) go to 



ip(x,t) — exp(— b c x) exp 



ifi 
2m 



b 2 c (t + t ) 



(33) 



which are formal solutions of the standard linear Schrodinger equation but are, evidently, physically unacceptable 
because they are not normalizable (and the wave function's square modulus |^| 2 itself diverges when xb c — > — oo). It 
is nevertheless interesting that the non-linearity associated with q > 1 not only regularizes (in the sense of leading to a 
finite norm N) the plane wave solutions of the Schrodinger equation (as stressed by NRT in [T]) but it also regularizes 
the exponential solutions (33) (at least for all times t > —to). 



E. (Non-)Preservation of the Norm 

It is known that, in general, time dependent solutions to the NRT nonlinear Schrodinger equation do not preserve 
the norm [2] . The g-plane wave solutions studied in [2] constitute a remarkable exception: they do preserve the norm. 
Up to now the g-plane wave solutions where the only known time dependent solutions to the NRT equation. This 
means that no norm non-preserving explicit solution was known before our present work. Therefore, it is of some 
interest to explore whether the solutions to the NRT equation investigated here preserve the norm or not. The norm 
N of the g-Gaussian wave packet is given by, 



N 



dx\^\ 2 , 



(34) 



where 



|V| 2 = [1 - 2(1 - q)Vl(a(t)x 2 + b(t)x + c(t)) + (1 - qf\a{t)x 2 + b(t)x + c(t)\ 2 



(35) 



When a^O the g-Gaussian wave packet is normalizable (that is, N < oo) provided that 1 < q < 5 and the polynomial 
P(z) = 1 — (1 — q)(a(t)z 2 + b(t)z + c(t)) doesn't have real roots. Notice that for a ^ the range of g-values admitting 
normalizable g-Gaussian wave functions of the form ^\ is larger than the range of g-values leading to normalizable 
g-plane wave functions (which is 1 < q < 3 [T]). 

Let us now consider a concrete instance of a solution to the NRT equation that does not preserve the norm. In the 
case of the solution (|3ll the norm is, 
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The norm is finite for t > —to and 1 < q < 3 and proportional to (f + to) 1-9 . Therefore, the norm is not conserved: 
it vanishes in the limit t — > +oo, it is a finite and monotonously decreasing function of time for all finite times t > t , 
and it diverges at t = —to- Consequently, (31) constitutes an explicit example of a norm non-preserving solution to 
the NRT equation. 

As a final comment, let us mention that the non conservation of the norm suggests that, when considering pairs of 
time dependent solutions tpi and tp2 to the NRT equation, the overlap J ip\^2dx is not conserved either. However, 
there may still be some conserved measure of the "distance" or "fidelity" between pairs of time dependent solutions. 
The search for such a measure constitutes a line of enquiry that may shed new light on the nature of the NRT 
dynamics. A possible direction to explore in this regard would be the one suggested by Yamano and Iguchi in 32], 
where non-Csiszar /-divergence measures where considered in connection with nonlinear Liouville-type equations. 
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III. WAVE PACKET SOLUTIONS IN THE PRESENCE OF AN HARMONIC CONFINING 

POTENTIAL 

Now we are going to use the g-Gaussian ansatz ([7]) to investigate time dependent solutions of nonlinear Schrodinger 
equation given by 



dt 2-q2m dx 2 K ' 



(37) 



where V(x) = ^Kx 2 . In the limit q — > 1 the above equation reduces to the time dependent Schrodinger equation of 
the quantum harmonic oscillator. Inserting the ansatz ([7| in the right hand side of equation (37) we get 



1 K Z d'V' q rr, Ma h 
— q 2^^x^ + V{x) ^ = 2^ 

~ 2m 



2o(t)(l - (1 - q)(a(t)x 2 + b(t)x + c(t))) - (2a(t)x + 6(i)) 2 J ip q + -Kx 2 ^ q 



K - 2(3 - q)a(ty x z - 2(3 - q)a(t)b(t)x + 2a(t) - 2(1 - q)a{t)c{t) - b(ty 



V> 9 



Comparing Eq.([8]) with Q we obtain the following set of coupled differential equations for the evolution of the 
parameters appearing in the ansatz 



ia(t) 
ib(t) 
ic(t) 



-(3- 9 )a(t) 2 -- 
m la 

-(3 - g)o(t)ft(*) 
m 

U(l-q)a(t)c(t)-a(t) + b{ f- 



(38) 



Integrating the first of these equation one obtains, 



1 f da{t) 
where 5 is an integration constant. 
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(39) 



A. Quasi-Stationary Solutions 



A particularly interesting solution of the set of equations ( 38 1 is given by, 



a = a c 
b = 0, 
c 



1 / mK 
nV 2(3 -q) ' 



1-9 



1 — exp — q) 



ha r t 



leading, in turn, to the following solution for the non-linear Schrodinger equation, 



(40) 



exp ( - g) 7 ^ ) - (1 - q)a c x 



(41) 



It is interesting to consider now the q — > 1 limit of the above solution, 



lim ib 

q-H T 



exp 
exp 
exp 



ha c t , 
-i hm 



exp ( — a c x 



1 - (1 - g)a c exp ( i(l - g) 



ha r t 



in 



exp (" 



raw 

2/T 



(42) 



which is the (unnormalized) wave function associated with the ground state of a standard harmonic oscillator with 
natural frequency uj = \ ^ and zero point energy E = \hio. 



The norm of the solution (41 1 is finite for 1 < q < 5 and — t c < t < t c with t c 



sigularities at t — ±t c . The time derivative of the norm of the solution (41 1 is 



(q—l)ha c 



and it has finite-time 



— = 2(1- q) sin 

at m 



(1-9) 



ha r t 



1 — 2a c (l — q) cos 



(1-9) 



ha r t 



dx. (43) 



The expression between square brackets appearing in the integrand in the right hand side of the above equation is in 
general larger than zero and, consequently, the time derivative of the norm is different from zero. Therefore, the time 
dependent solution (41) constitutes another explicit example of a solution that does not preserve the norm. In the 



limit q — > 1 we get iCr^ 00 an d, of course, dN/dt 0. 



IV. GENERALIZATIONS OF THE NRT APPROACH AND UNIQUE FEATURES OF THE NRT 

EQUATION 

The existence of (/-plane wave solutions consistent with the de Broglie relations connecting frequency and wave 
number respectively with energy and momentum was one of the features of NRT equation discussed in [T]. It is 
worth noticing that there are other possible nonlinear Schrodinger-like equations exhibiting similar properties. That 
is, the NRT approach can be substantially generalized. Let us consider a pair of one-variable functions L(u) and F(u) 
satisfying the functional relation, 



It can then be verified after some algebra that the (in general non linear) Schrodinger -like equation 



(44) 



dt 



$0 



2m dx 2 



L 



admits exact time dependent plane wave-like solutions of the form, 



$(M) 
$ 



(45) 



$(x,t) = <P Q F[i(kx - wt)], 



(46) 



with hw = h 2 k 2 /2m. Indeed, inserting the ansatz (461 into the right and left hand sides (45) and setting u = i(kx-wt) 
we get, 



if i 



dt 



2m dx 2 



L 



dF\ fdu 
~du) \dt 

h^d?_ 

2m du 2 



= hw 



dF 



[L(t 



du ' 
du s2 
dx 



h 2 k 2 d 2 

2m du 2 



[L{F{u))\ 



(47) 



where, in the last equation, the fact that d 2 u/dx 2 = was used. It is plain that the functional relation (44) implies 
(45), provided that hw = h 2 k 2 /2m. This in turn leads, via the de Broglie relations, to the kinetic energy-momentum 
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relation E = p 2 /2m. Similarly to what happens within the NRT scenario, the solutions (46 ) to the non linear evolution 



equation (45) propagate without changing shape and with constant velocity v = w/k. Furthermore, as we have seen, 



they comply with a relation between w and k that is consistent with the de Broglie connection between frequency, 
wave number, energy, and momentum. 



One procedure to generate pairs of functions satisfying (44 1 is the following. One starts with a function G{u) such 
that its derivative dG/du admits an inverse. Then we define, 



F(u) = 



dG 

du ' 



L{u) = G F^(u) , 

where F^~ x \u) is the inverse function of F(u), satisfing f (F(u)) = u. 
by (|48| comply with the required relation ( 44 ) . 



In the case of the NRT equation we have, 



(48) 



It can be verified that the functions defined 



G(u) 
F(u) 
L(u) 



' [1 + (I - q)u} 2 ^ . 



2-q 
[l + (l-q)u]T^ 



,2-9 



2-q 



(49) 



An example (different fro m t he NRT one) of a nonlinear Schrbdinger -like equation of the form ( 45 ) admitting the 
plane wave-like solutions (46) is given by F(u) = sinh(u) and L{u) = \/l + u 2 . This case corresponds to the nonlinear 
evolution equation, 



lh dt 



2m dx 2 



1 21 



$0 



(50) 



having plane wave-like solutions 



$(x,i) = $o sinh[i(fcc — wt)], 



(51) 



with hw — h 2 k 2 /2m. 

As mentioned in the Introduction, the NRT equation was inspired by the nonextensive generalized thcrmostatistical 
formalism based on the constrained optimization of Tsallis' entropy 1J. Indeed, there are formal connections between 
the NRT scenario and the alluded thermostatistical formalism. An intriguing question, which is beyond the scope 
of the present work but certainly deserves to be explored, is the possible existence of connections between (some of) 
the nonlinear Schrbdinger equations (45), on the one hand, and other formalisms based on non-standard entropic 



functionals different from the Tsallis one [U [33] , on the other hand 



We have seen that the NRT equation is not the only nonlinear equation of the form ( 45 1 admitting plane wave-like 



solutions consistent with the de Broglie relations. It is then natural to ask which of the equations (45) also admit 
solutions of the form 



${x,t) = <$> F[d 2 {t)x 2 + dx(t)x + d (t)], 



(52) 



with d\ (t) , di (t) , do (t) appropriate time dependent coefficients. The solutions (52) would constitute generalizations of 
the (/-Gaussian solutions to the NRT equation previously discussed in the present work. It turns out that the NRT 
equation is the only member of the family ( [45] ) admitting both solutions of the form ( 46 1 and of the form ( 52 ) . If one 

^- given by the ansatz 



inserts in the left and rigth hand sides of the evolution equation ( 45 ) the expression for tp 



ansatz ( 52 ) one gets 



ih 



<9V 

di 



ih(d 2 x 2 + dix + d )F'(u), 



(53) 
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and, 



2m 9a; 2 



L(F(u))} = - — \(2d 2 x + d 1 ) 
J Zm I 



du 2 



L(F(u)) 



2d Q 



du 



L(F(u))]}, 



(54) 



where u = d2(t)x 2 + d\(t)x + do(t). Then, to have solutions of the form (46) (plane wave-like), equation (44) must 
hold. On the other hand, to have solutions of the form (52) one further condition is required, 



du 



L(F{u)) = ( ri +r 2 u)F'(u), 



(55) 



with r% and r 2 appropriate constants. It follows from (55) that 



L"(F{u))F\u) = r 2 . 



(56) 



Combining now equations ( 44 ) , ( 55 ) and ( 56 1 one obtains 



which leads to, 



r 2 F'(u) + (n + r 2 u)F"(u) = F'(u), 



F'(u) = F (n +r 2 u)r?. 



where Fq is an integration constant. Finally, we have, 



(57) 



(58) 



F(u) = F a ( ri +r 2 u)-2 +F X 



(59) 



involving one more integration constant -Fi. Making now the identification r 2 = l — q one sees that (up to multiplicative 
and additive constants) the form of the solution (59) coincides with (/-Gaussian wave packet. 



CONCLUSIONS 



We have obtained a new family of exact, analytical time dependent wave packet solutions to the nonlinear 
Schrodinger equation recently proposed by Nobre, Rego-Montero and Tsallis [TJ Our solutions have the form 
of a ^-exponential evaluated upon a quadratic function of the spatial coordinate x with time dependent coefficients. 
Therefore, these solutions have a q-Gaussian form. They extend and generalize the previously known solutions to 
the NRT equation. The solutions investigated here by us correspond, in the limit q — > 1 of the parameter q, to the 
Gaussian wave packet solutions to the standard linear Schrodinger equation. Our present wave packet solutions admit 
as a special particular case the g-plane wave solutions studied in pQ . In the present work we also discuss other families 
of nonlinear Schrodingcr-like equations, besides the NRT one, leading to a dynamics compatible with the dc Broglie 
relations. In this regard, it is remarkable that the existence of the Gaussian-like time dependent solutions investigated 
in this work is a unique feature of the NRT equation not shared by the abovementioned, more general, families of 
nonlinear evolution equations. 

We also obtained q-Gaussian wave packet solutions for the case of a harmonic potential. In the limit q — >• 1 these 
latter solutions reduce to Schrodinger's celebrated Gaussian wave packet solutions to the harmonic oscillator. As 
a particular case of the time dependent q-Gaussian wave packets associated with the harmonic potential we found 
a quasi-stationary solution yielding in the q — > 1 limit the wave function corresponding to the ground state of the 
quantum harmonic oscillator. 
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